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I.  INTRODUCTION. 


For  many  years  it  has  been  observed  that  liquid-filled  projectiles  have  a 
proclivity  for  unusual  flight  behavior,  often  being  unstable  even  though  the 
same  projectile  with  a  solid  payload  is  stable.  The  spin  of  the  projectile 
imparts  rotation  to  the  liquid  which  affects  the  dynamic  behavior  of  the 
shell.  Typically,  the  fluid  is  contained  in  a  right-circular  cylinder.  In 
the  problem  considered  here,  the  fluid  fills  the  cylinder.  The  coordinate 
system  and  notation  are  shown  in  Figure  1, 

The  fluid  motions  can  be  separated  into  two  classes:  solid-body  rotation 
and  spin-up.  The  meaning  of  the  former  is  obvious.  The  latter  referb  to  the 
transient  state  produced  by  changing  the  fluid  motion  from  solid  rotation  at 
one  spin  (which  may  be  zero)  to  that  at  a  greater  spin.  The  physics  of  spin- 
up  from  rest,  the  problem  to  be  considered  here,  is  different  in  some  impor¬ 
tant  respects  from  that  of  spin-up  from  a  state  of  solid-body  rotation.1’2 
In  particular,  spin-up  from  rest  is  inherently  a  nonlinear  problem,  whereas 
spin-up  between  two  finite  states  of  solid-body  rotation  can  be  treated  as  a 
linear  problem  if  the  change  in  rotation  is  small.  Spin-down  is,  in  general, 
an  unstable  flow  and  often  becomes  turbulent.  Relatively  little  work  has  been 
done  on  spin-down;  see  Reference  3  for  a  discussion  of  this  problem  and  some 
pertinent  references  on  experimental  work. 

The  basic  reference  on  spin-up  from  rest  is  Wedemeyer's  paper?4  our  work 
is  an  extension  of  his.  Wedemeyer's  model  was  based  on  an  order-of-magnitude 
analysis  rather  than  on  formal  matched  asymptotic  expansions.  A  more  formal 
treatment  was  given  by  Greenspan1  for  a  two-term  expansion.  In  both  cases  a 
phenomenological  approach  was  required  to  complete  the  theory. 

Wedemeyer's  model  yielded  a  nonlinear  partial  differential  equation  of 
the  diffusion  type  for  the  azimuthal  velocity,  V,  as  a  function  of  t^rne  and 
radial  coordinate  but  not  of  axial  coordinate.  Without  the  diffusion  terms 
the  solution  for  V  is  elementary;  Wedemeyer  omitted  them  in  his  approximate 
analysis.  The  main  objective  of  this  report  is  to  present  the  numerical 


1.  H.P.  Greenspan,  The  Theory  of  Rotating  Fluids,  Cambridge  University  Press, 
London  and  New  York,  1&6&. 

2.  E.R,  Benton  and  A.  Clark,  Jr.,  "Spin-Up,"  article  in  Annual  Review  of 
Fluid  Mechanics,  Vol.  6,  Annual  Reviews,  Inc.,  Palo  Alto  .''California , 

1 97T. 

3.  G.P.  Neitzel  and  Stephen  H.  Davis,  '  ergy  Stability  Theory  of  Decel¬ 
erating  Swirl  Flows,"  The  Physics  of  Fluids,  Vol.  23,  No.  3,  March  1980, 
pp.  432-437. 

4.  E.H.  Wedemeyer,  "The  Unsteady  Flow  Within  a  Spinning  Cylinder,"  BRL  Report 

No,  1225,  October  1963  (AD  431846).  Also  Journal  of  Fluid  Mechanics,  Vol . 
20,  Part  3,  1964,  pp.  383-399.  ~~~~  “ 
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solution  of  Wedemeyer's  equation.  We  shall  state  the  necessary  equations  and 
boundary  conditions  and  shall  describe  the  numerical  procedure.  Also, we  shall 
present  results  together  with  some  discussion  of  the  accuracy  of  the  solu¬ 
tion.  Applications  to  projectile  problems  will  be  emphasized. 

One  of  the  main  applications  of  spin-up  theory  is  In  the  study  of  the 
eigenvalue  problem,  l.e.,  the  determination  of  frequencies  and  decay  rates  of 
the  waves  In  a  rotating  fluid.5  Spin-up  theory  Is  also  required  to  solve  the 
forced  oscillation,  or  moment,  problem.  Previously  thi.  theory  was  used  to 
study  the  spin  decay  of  a  liquid-filled  projectile6  after  ejection  from  the 
gun. 


Some  appreciation  of  spin-up  effects  on  a  projectile  flight  can  be  gained 
by  considering  orders  of  magnitude.  Let  a  and  c  be  the  radius  and  half-height 
of  the  cylinder,  .0  the  spin  of  the  projectile  in  rad/sec  (say  at  the  muzzle), 
v  the  kinematic  viscosity  of  the  liquid,  and 


Re  «  ft  a2/v,  E  =  v/fi  c2, 


(1.1) 


the  Reynolds  number  and  Ekman  number.  Note  that  E  is  sometimes  defined  using 
the  height,  2c,  rather  than  c,  and  that  for  historical  reasons  we  shall  use  Re 
in  this  report  rather  than  the  more  conventional  E.  If  the  Ekman  layers, 
i.e.,  endwall  boundary  layers,  are  laminar,  the  characteristic  time  for  spin-up 
is 


ts  *  (2c/a)  R e1/2/n 
=  2/E1/2n  (sec). 


(1.2) 


This  expression  is  derivable  from  linear  spin-up  theory1’2  and  from  the 
Wedemeyer  solution  without  diffusion.  We  also  use  a  nondimensional  spin-up 
time  tz  =  n  t$.  For  Re  >  105,  approximately,  the  Ekman  layer  may  be  turbulent, 

in  which  case  the  characteristic  spin-up  time  can  be  estimated  by 


E$t  =  (28.6  c/a)  Re1/5/n.  (1.3) 

This  result  can  be  obtained  from  the  Wedemeyer  solution  without  diffusion  for 
turbulent  Ekman  layers.  Basically,  these  are  time  scales  for  the  spin-up 
process  and  do  not  necessarily  give  a  measure  of  the  closeness  of  the  process 
to  solid-body  rotation  which  can  be  measured  in  several  ways,  e.g.. 


5.  R.  Sedney  and  N.  Gerber,  "Oscillations  of  a  Liquid  in  a  Rotating  Cylinder: 
Part  II.  Spin-Up BRL  Technical  Report  ARBRL-TR-02489 ,  May  2983 . 

6.  C.W.  Kitchens ,  Jr. ,  N.  Gerber ,  and  R.  Sedney,  "Spin-Decay  of  Liquid-Filled 
Projectiles Journal  of  Spacecraft  and  Rockets ,  Vol.  IS ,  No.  6,  November- 
Decemiber  1978,  pp.  348-354. 


0 


the  time  for  the  angular  momentum,  flow  rate  across  a  meridian  plane,  or  other 
property  to  attain  a  value  within  10%,  say,  of  the  solid-body  value.  These 
measures  require  detailed  calculations  and  so  are  not  very  convenient.  A  rule 
of  thumb  is  that  sol  id -body  rotation  is  reached  at  about  4  l  after  an  impul¬ 
sive  start  of  the  cylinder.  s 


For  projectile  flight  these  characteristic  spin-up  timos_should  be 
compared  with  tfl ,  the  time  of  flight  of  the  projectile.  If  tg  or 

<<  t^.j ,  spin-up  effects  can  be  disregarded  and  solid-body  rotation  can  be 

assumed.  If  tg  or  tgt  »  £^/10,  or  larger,  spin-up  effects  probably  have  to 

be  considered.  To  put  these  estimates  in  perspecti ve,  consider  the  parameters 
for  two  155mm  projectiles.  For  one  of  these.  Case  1,  c/a  =  3.120,  Re  = 


4  x  lO4,  a  =  754  rad/sec  giving  t  =  1.65  sec.  For  Case  2,  c/a  ■=  5.200,  Re  = 

2  x  106,  «  =  754  rad/sec  giving  t$t  =  3.52  see.  For  both,  t^  =  40  sec,  and  the 

above  criteria  indicate  that  spin-up  effects  may  be  important,  especially  for 
Case  2.  The  extent  to  which  the  projectile  motion  is,  in  fact,  affected  by 
the  internal  liquid  motion  requires  a  solution  to  the  forced  oscillation,  or 
moment,  problem.  Reference  7  treats  the  moment  problem  for  solid-body 
rotation  flow.  The  present  report  is  restricted  to  the  unperturbed  time- 
dependent  axisyminetric  fluid  motion. 


II.  THF.  SPIN-UP  MODEL 


A.  Time  Scales. 


Consider  the  motion  of  a  fluid  which  fills  a  cylinder,  initially  at  rest, 
wiiich  is  "rapidly"  brought  to  a  constant  angular  velocity.  The  fluid  motion 
is  axisyminetric  and  time  dep.-ndent.  The  conditions  for  which  the  Wedemeyer 
spin-up  model 4  is  valid  can  be  presented  in  terms  of  the  time  scales 
involved.  It  was  derived  on  the  basis  of  an  impulsively  started  cylindrical 
container,  a  theoretical  concept  which  will  be  discussed  below;  for  now  this 
time  scale  is  taken  to  be  zero.  Only  the  laminar  Ekman  layer  case  will  be 
discu  .ed.  There  are  three  time  scales  that  govern  the  spin-up  process.  In 
nondimensional  form  these  are 


-  2nft 

£s  =  2  E-1/Z  Si-1 
£„  *  (a/c)2  E"1  si"1. 


7.  N.  Gerber  and  R.  Sedney,  " Moment  on  a  Liquid-Filled  Spinning  and  Nutating 
Projectile:  Solid  Body  Rotation BRL  Teohnioal  Report  ARBRL-TR-02470 t 
February  1983  (AD  A12S332). 
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The  first  of  these  Is  the  time  for  the  container  to  rotate  through  one 
revolution.  It  can  be  shown  that  £  Is  the  time  scale  for  the  formation  of 

quasi-steady  Ekman  layers  on  the  endwalls,  which  layers  are  of  thickness 
\/n 

0(v/n)  .  Actually, the  Ekman  layers  are  quasi-steady  after  about  2  radians 

of  rotation,  i.e,,  q£  ■  2.  Wedemeyer  used  this  result;  a  more  complete  study 

of  Ekman  layer  formation  Is  given  by  Benton.8  The  thickness  of  the  Lkman 

1/2  -1/2 

layers  Is  assumed  to  be  small  compared  to  either  c  or  a  so  that  E  or  Re 
Is  required  to  be  small.  The  second  time  scale  Is  (1,2),  anticipated  to  be 

the  characteristic  spin-up  time;  spin-up  occurs  for  £/£g  ■  0(1).  The  time 

scale  £v  Is  the  time  for  velocity  gradients  or  vorticlty  to  diffuse  viscously 

across  the  radius  of  the  cylinder.  For  c/a  «  0(1),  t  Is  an  order  of 

magnitude  greater  than  £g  so  that  viscous  diffusion  plays  a  small  role  in 

spin-up.  For  an  Infinite  cylinder  It  Is  the  only  mechanism  for  spin-up;  for 

large  c/a,  £y  can  become  comparable  to  or  less  than  £,.. 


For  an  impulsive  start,  the  conditions  for  the  validity  of  the  spin-up 
model  are 


£ 

E 


ft 

-1/2 


«  c 


„  «  £  , 

S  V* 


>>  1  and 


c/a  b  0(1) . 


B,  Approximation  to  Impulsive  Start, 

The  degree  to  which  an  impulsive  start  can  be  approximated  in  an  experi¬ 
ment  depends  on  the  system  parameters.  cor  our  application  the  angular 
acceleration  is  large  during  the  time  when  the  projectile  is  in  the  gun;  the 
small  deceleration  in  flight,  the  spin  decay,6  is  neglected.  For  artillery 
projectiles  this  acceleration  time  is  typically  0.020  sec.  A  time  history  of 
projectile  spin,  n  (£),  is  shown  in  Figure  2.  To  quantify  the  departure  from 
an  ideal  Impulse  several  approaches  are  possible;  here  time  scales  are  used. 

For  a  finite  angular  acceleration  of  the  cylinder  let  the  spin  be 
?!  (£),  as  in  Figure  2,  e.g.  The  time  t  is  used  as  a  reference.  A  char¬ 
acteristic  time  £c  is  introduced,  the  time  when  the  acceleration  becomes  zero. 
The  acceleration  or  impulse  time  scale  is  taken  to  be 


8.  E.R.  Renton,  "On  the  Flow  Due  to  a  Rotating  Disk,"  Journal  of  Fluiu 
Meohanios .  Vol,  24,  Part  4,  1966,  pp.  781-800 » 
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which  was  used  by  Weidman9’ 10  ;  tj  =  0  for  a  true  impulse.  Tor  comparative 
purposes  ?<  (t)  is  taken  to  be  linear  for  0  <  t  <  t  and  p  (t)  =  s2  for 
t  >  E^.  In  most  spin-up  experiments  enough  information  to  define  \i  (t) 

more  precisely  is  not  available;  a  linear  approximation  to  d  (t)  in  Figure  2 
is  reasonable. 

In  Table  1  the  above  time  scales  are  given  for  three  experimental 
arrangements  that  have  been  used  in  spin-up  studies.  (In  some  other  experi¬ 
ments  insufficient  information  was  available  to  define  the  time  scales.)  The 
first  is  the  gun  tube,  using  the  parameters  of  Figure  2;  the  second  is  the 
spin  generator  from  which  Wedemeyer4  obtained  observations  of  secondary  flow 
during  spin-up;  the  third  is  the  most  nearly  impulsive  case  in  Weidman's 
experiments,  Figure  2c  of  Reference  10. 


TABLE  1.  TIME  SCALES  FOR  IMPULSIVE  START,  MS 


t. 

E 

E 

I 

ft 

c 

Gun  tube 

5.6 

10.8 

18 

Spin  generator 

26.6 

22.2 

200 

Weidman 

218 

57.6 

5,200 

To  approximate  an  impulsive  start  for  the 
must  be  satisfied:  E .  <<  tc  ai  a  condition  on 
fj  <<  to  insure  that  the  impulse  time  scale 


spin-up  problem,  two  criteria 
the  mechanical  system  and 
is  small  compared  to  the  time 


9,  p.p.  Weidman,  "On  the  Spin-Up  and  Spin-Down  of  a  Rotating  Fluid:  Part  1. 
Extending  the  Wedemeyer  Model."  Journal  of  Fluid  Meohanios,  Vol.  77,  Part 
4,  1976,'  pp .  685-708. 

10.  P.D.  Weidman,  "On  the  Spin-Up  and  Spin-Down  of  a  Rotating  Fluid:  Part  2. 
Measurements  and  Stability ."  Journal  of  Fluid  Meahanioe,  Vol.  77,  Part  4, 
1976,  pp.  709-735 . 


II 


of  formation  of  the  quasi -steady  Ekman  layers.  Table  1  shows  that  the  gun 
tube  gives  the  best  approximation  to  an  Impulsive  start,  but  only  marginally 
according  to  the  criteria.  The  spin  generator  does  not  give  a  Impulsive  start 
and  Weldman's  case  Is  far  from  It. 

Since  the  Wedemeyer  model  assumes  an  Impulsive  start,  deviations  from  It 
as  shown  In  Table  1  would  require  modification  of  the  model.  The  necessary- 
analysis  has  not  been  done  for  either  small  or  large  departures  from  an 
Impulsive  start.  In  References  9  and  11  the  same  type  of  modification  of  the 
Wedemeyer  model  was  made  to  account  for  a  nonlmpulslve  start.  The  value  of 
rt  was  replaced  by  the  variable  n  (£).  This  should  be  a  reasonable 
approximation  for  an  almost  Impulsive  start.  In  Reference  11  its  validity  was 
checked  by  comparison  with  solutions  to  the  Navler-Stokes  equations. 


C.  The  Physics  of  Spin-Up, 

In  time  of  O(^)  quasi-steady  Ekman  layers  are  formed  on  the  endwalls. 

For  small  time  and  near  the  center  of  the  endwalls  these  layers  are  essen¬ 
tially  the  same  as  the  boundary  layer  on  a  steady  rotating  disk,  the  von 
Karman  disk  problem.  The  suction  exerted  by  these  layers  draws  external  fluid 
Into  the  F.kman  layers  and  imparts  rotation  to  this  fluid.  With  no  pressure 
gradient  acting,  the  fluid  spirals  out  to  larger  radii,  where  the  Ekman  layers 
eject  fluid.  This  fluid  Is  now  rotating.  This  is  the  haslc  mechanism  for 
spin-up. 

0 *  Spin -Up  Th eory  . 

The  notation  here  will  be  the  same  as  in  Reference  5.  Lengths,  veloci¬ 
ties,  pressure,  and  time  are  made  nondlmensional  by  a,  art,  prt2a2,  and 

rt”\  respectively,  where  p  is  the  liquid  density  and  rt  is  the  constant  spin 
rate  of  the  cylinder..  In  the  inertial  frame  cylindrical  coordinates  r,  e,  z 
are  used,  with  the  origin  of  z  at  the  center  of  the  cylinder,  and  velocities 
U,  V,  W,  respectively.  Dimensionless  time  is  t.  Derivatives  are  indicated  by 
subscripts.  Wedemeyer4  showed  that  the  flow  can  be  divided  into  two  regions: 
the  quasi-steady  Ekman  layers  at  the  endwalls  and  the  rest  of  the  flow,  called 
the  core  flow.  Wedemeyer  did  not  point  out  that  a  boundary  layer-,  l.e., 
Stewartson  layer,  must  be  inserted  at  the  cylinder  wall.  Starting  with  the 
Navler-Stokes  equations  for  axi symmetric  flow. 


U*  +  U*  U*  + 
t  r 


W *  U*  -  V*2/r  =  -P*  +  Re' 
z  r 


( V2IJ*  -  U*/r2) 


(2.1a) 


n. 


C. IV.  Kitchens ,  Jr.,  and  N.  Gerber,  "Prediction  of  Spin-Decay  of  Liquid- 
Filled  Projectiles,"  BRL  Report  1996,  July  1977  (AD  A043275)  ■ 
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V*  +  U*  V*  +  W*  V*  +  U*  V*/r  =  Re'1  (v2V*  -  V*/r2)  (2.1b) 

W*  +  U*  W*  +  W*  W*  =  -P*  +  Re'1  V2  W*  (2.1c) 

r  z  z 

(rU*)r  +  rW*  =  0  (2. Id) 

V2()  =  ()rp  +  (1/r)  ()r  +  ()2Z,  (2.1e) 


where  the  asterisk  indicates  the  exact  solution,  he  used  order-of-magnitude 
arguments  as  indicated  in  Sections  !IA  and  IIC  to  simplify  them  in  the  core 
flow.  The  approximations  reduce  the  three  momentum  equations,  (2.1,  a,b,c), 
to 

Vt  +  U  (Vf  +  V/r)  >  Re-1  [Vpr  +  (V/r)f]  (2.2) 

and 

U2  =  Vz  =  Pz  =  0,  (2.3) 


where  the  asterisk  is  dropped  for  this  approximate  solution.  For  Re  ® 
Wedemeyer  proposed  neglecting  the  diffusion  terms  in  (2,2)  and  arrived  at 

Vw  +  Uw  (Vw  +  Vr)  =  °»  (2<4) 

t  r 

where  the  sub  w  indicates  this  approximation.  In  applying  his  model  Wedemeyer 
used  (2.4)  rather  than  (2.2).  Although  he  showed  the  crucial  importance  of 
the  Ekman  layers  to  the  spin-up  process,  his  model  did  not  require  a  solution 
for  the  flow  in  these  layers.  He  did  discuss  certain  properties  of  the 
solution  which  were  required  in  his  model  for  the  core  flow. 

A  mure  formal  approach  to  this  problem  was  given  by  Greenspan.1  In  his 
treatment,  lengths  are  made  dimensionless  by  2c  rather  than  a,  which  is  more 

natural  for  the  spin-up  process;  time  is  made  dimensionless  by  l  so  that  the 

1/2  b 

new  time  is  t'  =  t/t  =  t/[(2c/a)  Re  ].  An  expansion  in  the  small  parameter 
(l/ts)  =  ( a / 2 c )  Re"1/2  (2.5) 

is  assumed.  The  form  of  this  expansion 


U*  «  (1/t  )  Uj  +  .  .  . 

V*  -  v0  +  (i/ts)  vx  +  .  .  .« 

W*  =  (l/ts)  W1  +  .  .  . 

p*  =  P0  +  U/ts)  Pl  +  .  .  . 


(2.6) 
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follows  from  the  knowledge  that  In  the  core  flow  U*  and  W*  are  0(l/t  )  and  V* 

Is  0(1)  If  1/t  Is  small.  See  Reference  1  for  details.  The  Independent 

variables  are  r,  z,  t',  Substituting  the  expansion  Into  (2.1)  after  trans¬ 
forming  to  t'  yields 


V  +  ui  (V  v°/r)  30 


Vo2/r  -  Po 


=  vn  -  U.  =  0. 
°z  °z  *z 


(2.7) 


These  are  the  same  as  (2.3)  and  (2.4)  if  we  set  VQ  =  Vw,  Uw  =  (l/t$)  U^,  and 

t'  =  t/t  .  Although  the  formalism  of  matched  asymptotic  expansions  was  not 

used,  the  first  two  terms  of  the  outer  expansion  are  apparently  given  by 
(2.6).  The  Ekman  layer  solution  would  be  obtained  from  the  inner  expansion, 
which  was  not  considered  in  Reference  1. 


To  solve  (2.2),  (2.4),  or  (2.7)  a  relationship  between  U  and  V  is  neces¬ 
sary.  This  is  called  the  Ekman  compatibility  condition,  abbreviated  EC, 
because  the  Ekman  layer  suction  must  be  made  compatible  with  the  core  flow. 
Wedemeyer  used  the  facts  that  the  Ekman  layers  are  quasi -steady  after  one 
revolution  and  that  the  radial  mass  flux  in  the  core  flow  must  be  balanced  by 
that  in  the  Ekman  layers  to  obtain  some  conditions  on  the  U,  V  relationship. 

At  this  point  he  was  forced  to  take  a  phenomenological  approach.  The  rela¬ 
tionship  is  known  at  t  +  0  and  t  -*■  «  and  he  proposed  a  linear  interpolation 
between  them  to  obtain  an  approximate  relationship  for  any  t.  He  tested  this 
idea  in  some  other  problems  where  the  solution  was  known  and  decided  it  was 
satisfactory.  Some  confusion  has  appeared  in  the  literature  because  this  step 
was  misinterpreted;  further  discussion  of  this  is  given  below.  His  result  is 


U  «  k£  (V  -  r) 
k^  =  K  (a/c)  Re'1/2 


2</ts 


(2.8) 


for  laminar  Ekman  layers.  Wedemeyer  proposed  <  =  0.443  but  Greenspan  sug¬ 
gested  <  =  0.5;  the  latter  often  gives  results  in  better  agreement  with 
numerical  solutions  to  the  Navier-Stokes  equations.  For  turbulent  Ekman 


1  ayers 

U 


k 


t 


-kt  (r-V)8/5 
.035  (a/c)  Re"1/5 


(2.9) 


where  t  t  is  the  nondimensional ,  turbulent  spin-up  time.  The  core  flow  is 

assumed  to  be  laminar  so  that  turbulent  stresses  are  not  introduced  ,i n  the 
right-hand  side  of  (2.2). 
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Using  (2.8),  (2.4)  can  bo  solved  explicitly  for  V(r,t)  with  V(r,0)  -  0 
and  V(l,t)  *  1: 


2k  t  2k  t  -k  t 

V  =  (re  -  l/r)/(e  -  1)  for  r  >  e 

W 


0 


-V 

for  r  <  e  , 


(2.10) 


'k  £ 

so  that  r  =  e  separates  rotating  and  nonrotating  fluid;  it  is  called  the 
front.  Although  Vw  is  continuous  there,  a  discontinuity  in  shear,  i.e.,  Vw  . 

exists  because  diffusion  is  neglected.  Using  (2.9)  a  numerical  integration  is 
necessary,  but  the  character  of  the  solution  is  the  same  as  when  (2.8)  is 
used.  The  radial  velocity  is  obtained  from  (2.8)  or  (2.9)  and  W  is  obtained 
from  the  continuity  equation 


W  =  -U/r)  (rU)p 

At  r  =  1,  W  *  0  so  that  a  Stewartson  layer  must  be  inserted  there,  as  men¬ 
tioned  before;  this  has  yet  to  be  done.  Also  W  *  0  at  the  endwalls  z  =  +  c/a 
and  W  is  discontinuous  at  the  front. 


Wedemeyer  pointed  out  that  the  shear  discontinuity  in  the  solution  to 
(2.4),  see  (2.10),  would  be  smoothed  out  if  the  diffusion  terms  were  retained, 
as  in  (2.2).  But  this  is  a  nonlinear  second  order  equation  and  must  be 
integrated  by  finite  difference  methods.  This  is  discussed  in  Section  III. 

For  laminar  Ekman  layers,  the  spin-up  velocity  profile  is 

V  «  f  (r,  ktt,  kt  Re), 


where  k^t  =  t'  if  <  =  0.5;  for  the  turbulent  case  k^  is  replaced  by  k^.  These 

functional  forms  follow  directly  from  (2.2)  and  (2.8)  or  (2.9)  and  represent 
the  most  efficient  way  to  display  the  results. 


E.  Discussion  of  Ekman  Compatibility  Condition. 

The  original  forms  of  the  EC  given  by  Wedemeyer,  for  laminar  and 
turbulent  Ekman  layers,  were  discussed  in  Section  I  ID .  It  is  important  to 
realize  that  an  exact  EC  does  not  exist;  it  is  necessarily  approximate.  This 
can  be  deduced  theoretically  and  is  shown  from  finite  difference  solutions  of 
the  spin-up  problem  using  the  Navier-Stokes  equations.12  Three  attempts  to 
obtain  EC  which  give  improved  approximations  to  V  will  be  described. 


12.  C ,W .  Kitchens ,  Jr.,  "Ekman  Compatibility  Conditions  in  Wedemeyer  Spin-Up 
Model,"  The  Physics  of  Fluids,  Vol.  23,  No. 5,  May  1980,  pp.  1082-1064 . 
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From  his  finite  difference  solutions  Kitchens12  obtained  the  radial  mass 
flow  In  the  Ekman  layer  as  a  function  of  V/r  for  several  r  and  t.  The  mass 
flow  can  be  used  to  form  the  EC.  A  unique  relation  was  not  found,  rather  a 
band  of  poIntSc  Based  on  these  numerical  data.  Kitchens  deduced  an  approxi¬ 
mate,  monotonlc,  nonlinear  EC.  For  0  <  V/r  <  0.75,  U  Is  a  cubic  In  V/r  and 
for  0.75  <  V/r  <  1,  it  Is  a  linear  function.  Limited  tests  of  this  EC  show 
only  small  improvement  in  V  compared  to  using  (2.8). 

The  following  EC  is  proposed  as  an  extension  of  (2.8): 

U  =  -[r/(Re1/2  c/a)]  [.500  (1-V/r)  -  .057  (1-V/r)3]  (2.11) 

This  relation  satisfies  (9)  in  Reference  4  which  holds  for  V/r  =  0  at  t  =  0 
and  the  condition  U  (V/r  *  1)  =  0,  for  t  ■*■  »,  the  two  conditions  imposed  by 
Wedemeyer  in  arriving  at  his  linear  expression.  However,  it  also  satisfies 
(14)  in  Reference  4,  in  slope  and  curvature,  when  the  core  flow  is  almost 
solid-body  rotation.  Only  one  test  of  this  EC  has  been  made;  It  showed  that  V 
was  significantly  closer  to  the  results  from  Navier-Stokes  calculations  than 
those  using  other  EC's.  This  EC  is  the  only  one  which  is  a  rational  extension 
of  Wedemeyer' s  linear  expression. 

The  last  EC  to  be  discussed  is  of  a  somewhat  different  nature,  and  its 
origin  is  not  clear.  It  utilizes  the  solutions  of  Rogers  and  Lance13  to  the 
boundary  layer  flow  over  an  infinite  disk  rotating  with  angular  velocity  .0  and 
outer  flow  in  solid-body  rotation  with  angular  velocity  w  .  The  main  result 
used  from  this  paper  is  the  computed  radial  mass  flux  in  the  boundary  layer/ 

r  Re-1^2  =  m  vs  u>/2.  It  Is  plotted  in  Figure  3,  and  we  shall  call  this  the  RL 
curve.  Using  arguments  somewhat  different  from  Wedemeyer's,  Greenspan1  de¬ 
rived  the  EC  (2.8)  using  a  linear  approximation  to  the  RL  curve.  Possibly 
because  of  a  suggestion  of  Greenspan,1  page  169,  several  investigators  have 
attempted  to  improve  (2.8)  by  using  an  analytical  fit  to  the  RL  curve;  these 
include  Goller  and  Ranov,14  Venezian15  Watkins  and  Hussey,16  and  Weidman.9 
This  EC  is  determined  from  Figure  3  .sing  the  alternate  labels  on  the  scales, 

U  (c/a)  Re  /r  vs  V/r.  Weidman9  uses  a  7th  order  polynominal  to  approximate 
the  RL  curve. 


13.  Mr  H.  Rogers  and  G.N.  Lance,  "The  Rotationally  Symmetric  Flow  of  a  Viscous 
Fluid  in  the  Presence  of  an  Infinite  Rotating  Disk ,"  Journal  of  Fluid 
Mechanics ,  Vol.  7,  Part  2,  1960 ,  pp.  617-631 . 

14.  H.  Goller  and  J.  Ranov ,  "Unsteady  Rotating  Flow  in  a  Cylinder  with  a  Free 
Surface Journal  of  Basic  Engineering ,  TRANS.  ASMS,  Vol.  90 ,  Series  D, 
December  1968,  pp.  445-454. ~ 

15.  G.  Venezian ,  Topics  in  Ocean  Engineering ,  Vol.  2,  pp.  87-96,  Gulf 
Publishing  Co. ,  Houston,  Texas,  1970. 

16.  W.B.  Watkins  and  R.  G.  Hussey,  "Spin-Up  from  Rest:  Limitations  of  the 
Wedemeyer  Model,"  The  Physics  of  Fluids,  Vol.  16,  No.  9,  September  1973, 
pp.  1530-1531. 
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Note  that  Figure  3  appears  in  Wedemeyer' s  paper'1  hut  rot  tor  the  pu rnoso 
ot  determi ring  an  EC.  After  presenting  his  arguments  fur  (2.8),  he  tested  it. 
against  two  other  similar  problems  with  known  solutions.  One  of  the"."  was  the 
solutions  of  Rogers  and  Lance,  l 3  There  has  been  some  confusion  on  this 
matter  because  it  was  assumed  that  Wedemeyer  obtained  (2.8)  as  a  linen' 
approximation  to  the  RL  curve. 

As  pointed  out  by  l.c-idman,  one  of  the  assumptions  needed  to  apply  the 
Rogers  and  Lance  results  to  derive  an  EC  is  that  the  fluid  in  the  core  is 
"locally  in  solid-body  rotation."  It  appears  tnat  this  assumption  is  rather 
unrealistic.  Some  arguments  can  be  constructed  to  show  that  this  EC  is 
implausible;  but  it  cannot  be  proved  that  it  is  incorrect,  because  an  exact  FC 
does  not  exist.  Some  qualitative  arguments  were  given  by  Renton;17  later,  he 
provided  some  arguments  of  a  more  quantitative  nature18  which  were  based  on 
the  computations  described  in  this  report.  One  can  discuss  the  consequences 
of  using  the  EC  obtained  from  the  RL  curve  for  either  (2.2),  the  spin-up 
equation  with  diffusion,  or  (2.4),  without  diffusion.  Some  examples  of 
results  from  the  numerical  solution  of  (2.2),  using  this  EC,  will  he  given 
later.  The  solution  of  (2.4)  with  the  EC  using  the  RL  curve  has  some  unusual 
properties  which  actually  provide  the  most  compelling  reason  for  rejecting 
this  EC.  Since  (2.4)  is  a  first  order,  quasi-linear  partial  differential 
equation,  it  can  be  solved  by  quadrature;  this  is  true  if  (2.9)  is  used, 
e.g..  For  this  purpose  the  method  of  characteristics  is  used,  there  being  a 
one  parameter  family  of  characteristics.  Using  the  EC  from  the  RL  curve, 
Weidinan9  showed  that  the  characteristics  intersect  near  the  front,  implying 
that  the  solution  to  (2.4)  is  douhle-valued:  for  a  fixed  t,  V  is  double¬ 
valued  over  a  range  of  r,  in  the  neighborhood  of  the  r  for  which  V  =  0. 

He  showed  that  this  result  was  directly  attributable  to  the  fact  that  the  RL. 
curve  is  not  monotonic.  Clearly  this  is  an  unacceptable  result.  Weidinan  says 
that  there  must  be  a  discontinuity  in  V  to  resolve  this  dilemma.  This  remains 
a  conjecture  since  he  did  not  determine  the  position  and  strength  of  the 
discontinuity. 

Using  the  EC  from  the  RL  curve  in  the  spin-up  equation  with  diffusion, 
(2.2),  the  numerical  solution  is  not  double-valued  or  discontinuous.  One  such 
result  will  be  shown  later.  It  will  be  seen  that  when  the  limit  with  no 
diffusion  is  approached,  Re  -*•,  »,  there  is  a  tendency  for  the  V  vs  r  curve 

to  develop  a  vertical  tangent.  This  might  he  an  indication  of  an  approach  to 
a  discontinuity.  Both  a  vertical  tangent  in  the  theory  with  diffusion  or  a 
discontinuity  without  diffusion  are  physically  unacceptable  and  violate  the 
assumptions  of  the  model.  Since  the  double-valued/discontinuous  nature  of  the 
solution  is  caused  by  the  nonmonotonic  RL  curve,  it  seems  evident  that  the  EC 
based  on  it  should  be  discarded.  Any  suggestion  that  these  anomalies  indicate 
a  breakdown  of  the  Wedemeyer  model  is  unjustified,  of  course.  The  numerical 
solutions  of  (2.2),  rather  than  (2.4),  using  this  EC  does  give  reasonahle- 


17.  E.R.  Benton,  Memo  to  COTR ,  Taek  Order  74-461,  BRL,  January  1975. 

18.  E.R.  Benton,  "V ortioity  Dynamics  in  Spin-Up  from  Rest,"  The  Physics  of 
Fluids,  Vo l.  22,  No.  7,  July  1979,  pp.  1250-1251. 
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looking  V(r,t)  as  long  as  the  limit  with  no  diffusion  Is  not  approached,  1 .e. L 
If  k  Re  Is  not  large.  For  reasons  given  above  this  EC  Is  not  used  except  for 

showing  some  comparative  results. 


III.  METHOD  OF  SOLUTION 
A.  Finite  Difference  Method. 


The  initial  and  boundary  conditions  for  (2.2)  are: 


V  *  0 

for  t  <  0, 

0  <  r  <  1 

(3.1) 

V  =  1 

for  r  =  1, 

t  >  0 

(3.2) 

V  -  0 

for  r  *  0, 

t  >  0 

(3.3) 

Since  (2.2)  is 

a  parabolic,  quasi 

-linear  equation, 

it  can  be  solved  by  march- 

ing  in  time;  i.e.,  V(r,t)  is  obtained  at  t  =  kAt,  where  k  =  1,  2,  .  .  . 
and  At  is  the  time  increment  used  in  approximating  V^,.  The  grid  used  in  the 

finite  difference  scheme  is  shown  In  Figure  4.  The  scheme  Is  an  implicit  one 
and  the  discretization  error  is  second  order  in  both  r  and  t,  since  central 
differences  are  used  in  approximating  all  derivatives. 

The  interval  0  <  r  <  1  is  divided  into  N  segments  of  equal  length 
Ar  =  l/N.  V  is  calculated  at  the  N-l  nodal  points  r^  =  j  Ar,  j  *  1, 

2,  .  .  .,  N-l.  Given  V  at  t  *  kAt,  for  all  j,  the  solution  at  (k+l)At  Is 
obtained.  At  that  time  level,  consider  a  typical  point  Q  in  Figure  4.  The 
finite  difference  approximation  to  (2.2)  is  centered  at  P.  Thus. 


Vq  -  VR  =  (1/2)  (Fr  +  Fg)  At  +  0  (At)3  (3.4) 

where 

F  =  Re'1  [Vrr  +  (V/r)r3  -  U  (Vp  +  V/r)  (3.5) 

The  following  difference  approximations  are  made  at  the  point  R  (*  jAr,  kAt): 

vr  =  (VB  -  VA)/( 2  Ar)  (3.6a) 

Vrr  =  <VB  "  2VR  +  VA)/(Ar)2  (3.6b) 

(V/r) p  -  [(V/r)B  -  (V/r)A]/(2  Ar)  (j  >  1)  (3.6o) 

(V/r)p  =  (1/2)  (VB  -  2VR)/(Ar)2  (j  =  1).  (3.6d) 


The  special  form  of  (3.6d)  is  required  because  (V/r)A  is  indeterminate  at 

r  =  Ar.  Identical  approximations  are  made  at  point  Q,  with  points  A  and  R 
replaced  by  points  S  and  T,  respectively,  in  (3.6). 
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In  (3,4)  arui  (3.5)  IJq,  related  to  Vg  by  the  EC,  requires  an  iterative 

process.  An  initial  guess  for  V(,  is  used  to  calculate  Ug.  The  calculated  Vq 

fran  (3.4)  is  then  used  to  obtain  a  new  IL,  etc.,  until  th  'sired  degree  of 

convergence  is  obtained  at  all  points  for  t  °  (k  +  1)  At.,  n'tial  guess 

for  Vq  is  obtained  by  extrapolating  from  and  V^.  Recause  of  the  implicit 

nature  of  the  finite  difference  scheme  the  values  of  V  at  all  nodal  points  at 
t  =  ( k + 1 )  At  must  be  obtained  simultaneously.  However,  the  matrix  of  the  set 
of  N  linear  algebraic  equations  is  tridiagonal  so  that  a  straightforward 
algorithm  can  be  used  to  solve  the  equations,  which 


aj  Vj-1,  k+1  +  bj  v.j,  k+1  +  Lj  vj+l,  k+1 
The  expressions  for  a-,  b-,  c.,  and  d-  are  given  in 

%}  J 

3.  Accuracy  of  Results. 

Convergence  criteria  were  set  to  insure  that  the  numerical  solution  was 
accurate  ai  at  least  three  decimal  places.  For  a  number  of  cases,  covering  a 
large  range  of  Re,  or  Re,  calculations  were  made  for  different  com¬ 
binations  of  at  and  At  to  determine  maximum  permissible  interval  sizes.  The 
At  was  kept  small  enough  so  that  no  more  than  four  iterations  were  required 
for  convergence  at  any  time  step.  The  Ar  =  1/N  required  to  compute  V  to  an 
accuracy  of,  say,  1%  or  less  depends  on  r  and  t  and  on  the  parameter  kfc  Re,  or 

kt  Re,  Some  representative  results  showing  the  accuracy  of  V  are  presented  in 

Tables  2  and  3  for  the  parameters  of  Cases  1  and  2,  respectively.  The  per¬ 
centage  error  is  largest  for  small  r;  V  at  r  =  0,2  is  chosen  for  Tables  2  and 
3.  For  k ^  Re  =  23  in  Table  2,  N  =  100  is  large  enough  for  an  accuracy  of  1% 

or  less  when  t  >  600.  For  k.  Re  -  736  in  Table  2,  N  >  200  is  required  for 

t  =  6,000.  When  this  solution  is  used  as  the  basic  flow  in  our  eigenvalue 

calculations,  the  vorticity  t  =  (rV)  /2r  is  also  required.  The  error  in  it  is 

usually  largest  near  r  =  1.  A  convenient  test  on  the  accuracy  of  c  is  avail¬ 
able  and  discussed  in  Section  III  C.  The  CPU  time  for  solving  (2.2)  varies 
with  the  parameter  but  is  usually  less  than  one  minute  on  the  VAX. 

Ultimately,  the  adequacy  of  the  numerical  solution,  and  of  the  model 
itself,  must  be  judged  by  comparisons  with  solutions  to  the  Navier-Stokcs 
equations.  Some  examples  of  these  comparisons  will  be  given  in  the  next 
section. 


have  the  form 
d i .  (3,7) 

the  Appendix. 


TABLE  2.  EFFECT  OF  N  'N  SPIN-UP  VELOCITY  PROFILE.  CASE  1,  k*  Re 


t  -  1,200, 

r  ■  0,20 

N_ 

V_ 

av/v  « 

[  V ( N) 

-  V(100)]/V(100) 

20 

.07558 

3.49% 

40 

,07355 

0.71% 

80 

.07309 

0,08% 

100 

.07303 

0.00% 

t  =  600, 

r  -  0.20 

_N 

V_ 

AV/V  = 

[  V  ( N) 

-  V(100)]/V(100) 

20 

.00240 

72.6% 

40 

.00160 

.5.2% 

80 

.00141 

1.4% 

100 

.00139 

0.0% 

TABLE  3.  EFFECT  OF  N  ON  SPIN-UP  VELOCITY  PROFILE,  CASE  2,  k«.  Re 


t  =  9,000, 

r  =  0.20 

N 

V_ 

A V/V  * 

[  V  ( N) 

-  V(400)]/V(400) 

50 

.08626 

1.90% 

100 

.08517 

0.63% 

200 

.08476 

0.14% 

400 

.08464 

0.00% 

t  *  6,000, 

r  ■  0.20 

1 

_V 

A  V/V  » 

[V(  N) 

-  V( 1000) 3/V( 1000) 

50 

.01748 

21.2% 

100 

.01534 

6.4% 

200 

.01461 

1.3% 

400 

.01440 

-0.14% 

800 

.01443 

0.07% 

1000 

.01442 

0.00% 

0 


_ 
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For  an  imoulsive  start,  conditions  (3.1)  and  (3.?)  require  a  discon¬ 
tinuity  in  V  at  r  =  1,  t  =  0.  If  the  finite  difference  scheme  just  described 
is  used  without  modification,  the  value  of  V ( 1 ,0 }  would  be  required  but  is  not 
available.  If  values  such  as  0,  1,  or  1/2  are  used,  a  significant  error  is 
made  in  a  small  neighborhood  of  r  s  1.  Because  of  the  diffusive  nature  of 
(2.2),  this  error  will  decrease  as  t  increases.  The  error  may  be  significant 
at  the  earliest  time  for  which  a.  solution  is  required.  To  reduce  this  error 
a  local  solution  is  needed  to  resolve  the  discontinuity. 

Introducing  new  coordinates 

R  =  (1  -  r)  t'“1/2,  T  =  t' , 

a  solution  for  small  t'  was  obtained  in  the  form 

V  =  Vo  (R)  +  T1/2  Vj  (R)  +  .  .  . 

The  first  term  is  given  by 

VQ  =  erfc  (R/2>:^2)  =  erfc  ([1-r]  [Re/t] ^^2/2 ) 

=  l/(kA  Re). 

For  r  near  unity,  this  is  the  same  as  the  solution  for  the  impulsively 
started,  infinite  plate  (the  Rayleigh  problem)  as  expected.  This  solution 
satisfies  (3.1)  and  (3.2)  but  not  (3.3).  It  is  used  to  provide  the  initial 
condition  which  is  applied  at  the  first  time  step  t  =  At.  Since  (3.3)  is  not 
satisfied,  ari  error  is  made  at  r  =  0.  A  limit  is  placed  on  this  error:  V 

(0,  At)  <  lCf6,  which  is  satisfied  if  (Re/At)1//c  >  7.  The  latter  inequality 
must  be  satisfied  for  the  calculation  to  proceed.  In  practice  it  has  never 
been  violated.  The  second  term,  7^,  has  not  been  needed  in  our  calculations. 

Using  this  initial  condition  we  obtain  V/r  vs  r  shown  in  Figure  5  for 
Re  =  4974,  c/a  =  3.30  and  t  -  20  and  50.  Also  shown  are  the  results  obtained 
if  V ( 1 ,0 )  =  1  is  assumed;  there  is  a  relatively  large  error  for  t  =  20.  The 
intersection  of  the  dashed  lines  is  a  consequence  of  the  error.  In  Figure  6 
the  nondimensional  vorticity  r,  -  (rV)  / 2r  vs  r  is  shown  for  t  -  100  and 
600.  The  two  sets  of  curves  are  againobtained  using  this  initial  condition 
and  V(1,0)  =  1.  For  r  >  0.95  and  t  =  100  the  errors  resulting  from  use  of 

V(1,0)  =  1  are  quite  large.  The  vorticity  plot  has  the  advantage  that  a 
check  is  readily  available.  It  can  be  shown  from  (2.2)  and  (2.3)  that 
C  (r  =  1)  =  0.  The  slope  of  the  solid  curve  at  r  =  1  for  t  =  100  is  slightly 

negative;  for  t  =  600  it  is  zero.  The  slopes  of  the  dashed  curves  at  r  =  1 
are  quite  different  from  zero. 


IV.  RESULTS 


The  main  features  of  the  theo.y  win  be  Illustrated  by  some  examples. 
Usually  the  parameters  are  chosen  to  be  of  Interest  In  liquid-filled  p-ojec- 
t 11 e  applications,  but  sometimes  they  are  chosen  to  make  a  particular  point. 

The  Wedemeyer  model  for  spin-up  simplifies  the  Navier-Stokes  equations  so 
that  a  more  tractable  mathematical  problem  is  obtained  while  retaining  the 
essential  physics  of  the  problem.  Since  It  is  now  feasible  to  solve  the 
Navier-Stokes  equations  by  finite  difference  techniques,  it  may  be  that  the 
model  is  not  needed;  it  Is  premature  to  draw  that  conclusion.  However,  the 
results  of  the  model  can  be  compared  with  accurate  finite  difference  solu¬ 
tions.19  One  such  comparison  is  shown  in  Figure  7.  Using  the  program  of 
Reference  19,  V  vs  r  was  computed  for  Re  =  103,  t  *  62.5  and  for  Re  =  1014, 
t  *  62.5  and  125,  both  for  c/a  *  1.  The  results  at  z  a  0  are  plotted.  The 
circled  points  are  solutions  of  (2.2)  and  (2, .8).  For  these  parameters  the 
differences  between  the  spin-up  model  and  the  Navier-Stokes  solutions  are  of 
the  order  of  a  few  percent,  except  possibly  for  very  small  V.  The  differences 
are  not  always  so  small;  v,*  have  no  general  rules  for  what  values  of  the 
parameters  give  such  small  differences. 

Consider  the  solution  to  the  spin-up  equation  for  the  parameters  of  the 
two  cases  of  155mm  projectiles,  given  in  the  Introduction.  For  Case  1,  with 
Re  =  4  x  104,  we  should  expect  laminar  Ekman  layers.  The  solution  is  shown  in 

Figure  8  for  three  values  of  t.  One  of  them  is  t  =  1245.  For  t  =  4800, 

approximately  4  t  ,  solid-body  rotation  is  achieved,  essentially.  For  Case  2, 

with  Re  =  2  x  106,  we  should  expect  turbulent  Ekman  layers.  The  solution  to 
(2.2)  with  (2.9)  gives  the  velocity  profile  shown  in  Figure  9  for  three 
values  of  t  including  t  t  *  2700.  For  t  =  10,000,  approximately  4  tgt,  solid- 

body  rotation  has  not  been  achieved  to  the  same  degree  as  in  Figure  8, 

The  Ekman  compatibility  conditions  (2.8)  and  (2.9)  were  proposed  by 
Wedemeyer  for  either  laminar  or  turbulent  Ekman  layers,  there  being  no  pro¬ 
vision  for  transition  from  one  to  the  other.  Quite  different  results  con  be 
obtained  using  (2.8)  or  (2.9).  This  is  illustrated  in  Figure  10  with  Re  = 

3  x  105,  a  value  at  which  either  condition  might  apply.  The  results  using 
(2.8)  and  (2.9)  are  shown  for  two  values  of  c/a  and  t.  The  difference  is 
large  for  c/a  =  1.0  but  small  for  c/a  =  0.05,  an  extremely  small  aspect  ratio. 

The  need  for  including  diffusion  in  the  spin-up  model  was  discussed  in 
Section  II.  The  differences  in  V,  with  and  without  diffusion,  are  illustrated 
in  Figure  11  for  Re  ®  6.08  x  105,  c/'a  =  2.679  and  t  =  800  and  3200.  The 
differences  can  be  smaller  or  larger,  depending  on  the  parameters;  generally 
they  are  smaller  at  larger  values  of  k  Re.  For  given  t,  the  largest 
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difference  occurs  at  the  front  and  this  increases  as  t  decreases.  Even  though 
Re  is  large,  laminar  i.kman  layers  are  assumed  in  this  example. 

A  similar  comparison  is  shown  in  Figure  12  for  turbulent  Lkman  layers  for 
Re  =  2  x  10b  and  c/a  =  5.2.  The  results  including  diffusion  are  the  same  as 
those  in  Figure  9  for  t  =  1J00  and  2700.  Without  diffusion,  V  depends  only  on 
r  and  k^t.  The  differences  are  smaller  than  for  the  laminar  case. 

The  effect  on  the  velocity  profile  of  using  the  CC  from  the  Rl.  curve 
(Figure  3)  is  shown  in  Figure  13.  A  rather  abrupt  change  in  curvature  is 
obtained  for  small  V;  recall  that  the  maximum  in  the  RL  curve  occurs  for  small 
V.  To  obtain  this  result  it  was  necessary  to  use  ar  =  1/200;  a  value 
of  Ar  *  1/50  did  not  show  the  abrupt  change  in  curvature.  The  purpose  of  this 
calculation  was  to  illustrate  the  approach  to  the  no-diffusion  limit  in  which 
double-valued,  or  possibly  discontinuous,  V  is  predicted  by  Weidman's 
theory.9  Thus  the  rather  large  value  Re  =  5000,  It  is  not  unreasonable  to 

conclude  that  the  curve  in  Figure  13  is  tending  toward  one  with  a  vertical 
tangent,  which  might  be  related  to  the  predicted  douhl e-val ued  or  discon¬ 
tinuous  V.  As  the  slope  of  the  tangent  i ncreases,  the  calculation  would  become 
increasingly  more  difficult.  For  comparison  the  V  calculated  from  (2.2)  with 
(2.8)  is  also  shown.  The  difference  between  the  two  is  about  0.1  at  r  1  0,4, 

which  is  quite  large;  for  smaller  k  q  Re  the  differences  are  smaller.  The 
results  obtained  by  using  (2.8)  show  no  unusual  behavior.  This  comparison 
gives  further  evidence  for  not  using  the  EC  from  the  RL  curve. 


V.  APPROXIMATIONS  FOR  DIFFUSION  EFFECTS 

Inclusion  of  diffusion  effects  in  the  solution  to  the  spin-up  equation 
was  amply  demonstrated  in  the  previous  section,  but  the  numerical  results  do 
not  provide,  much  understanding  of  the  physical  process.  Without  diffusion  a 
shear  discontinuity  exists  at  the  front  which  implies  that  viscous  diffusion, 
neglected  in  arriving  at  (2.4),  cannot  be  neglected  in  the  neighborhood  of  the 
front.  A  local  solution,  including  diffusion,  can  be  sought  to  examine  the 
structure  of  the  shear  layer.  This  was  done  by  Venezian,19  who  did  a  boundary 
layer  type  analysis  of  (2.2)  on  the  moving  front.  His  result  will  be  shown 
below  as  part  of  a  more  general  solution.  He  found  that  the  moving  shear 

1/4 

discontinuity  is  a  layer  of  thickness  Q(E  ).  His  interpretation  is  that  the 
1/4 

E  Stewartson  layer  that  exists  at  the  sidewall  in  linear  problems  breaks 
away  from  the  sidewall  for  this  strongly  nonlinear  case  and  propagates  into 
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the  interior.  The  E  '  Stewartson  layer  apparently  remains  attached;  it  is 
the  one  mentioned  just  before  (2.1a). 

Inclusion  of  diffusion  over  the  entire  radius  can  be  studied  using 
matched  asymptotic  expansions.  It  is  convenient  to  work  with  ‘•he  circu¬ 
lation  r  3  rV.  Using  <  =  0.5  and  introducing  the  time  t',  (2.2)  and  (2.8) 
yield 


rt.  -  (r  -  r/r)  "r  =  £  (rrr  -  rr/r) 


(5.1) 
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where 


e  a  (k£  Re) 


-1 


Coordinates  centered  on  the  front  are  Introduced 

2t ' 


x 

y 


/  2  2t 
(re 

2t ' 

e^  -  1. 


l)/(e 


1) 


After  transforming  (5.1)  to  the  (x,y)  coordinates,  it  is  solved  by  expansions 
in  e.  The  first  term  in  the  outer  solution  is  (2.10);  the  first  term  in  the 
inner  solution  1$  Venezlan's  solution.  Using  a  2-term  Inner  expansion  and  a 
4-term  outer  expansion,  the  composite  solution  is 


r  *  e1/2  2( 2/ -rr) 1/2  y”1/2  e“S  /8  [erfc  (S/81/2) H  ~ 1  +  e  H  (S)  (5.2) 

where 

S  =  x(y/e)1/2. 

The  first  term  is  essentially  the  same  as  Venezlan's  solution.  H  (S)  satis¬ 
fies  a  linear,  2nd  order,  ordinary  differential  equation.  It  was  most  con¬ 
venient  to  solve  that  equation  numerically,  although  asymptotic  forms  were 
also  derived.  The  solution  H  (S)  is  shown  in  Figure  14. 

In  Figure  15  V  vs  r  is  shown  for  e  =  1/50  at  t'  =  0.80.  The  solution 
(2.10)  with  the  shear  discontinuity  is  shown  and  the  solution  to  (2.2)  with 
(2.8).  Venezian's  result,  the  first  term  of  (5.2),  is  also  presented;  his 
solution  is  very  successful  in  correcting  (2.10)  in  the  neighborhood  of  the 
front  and  for  larger  radii.  As  r  decreases  from  the  radius  of  the  front,  his 
solution  deviates  increasingly  from  the  solution  to  (2,2)  and  has  a  minimum. 

Of  course  it  was  intended  to  apply  only  near  the  front.  Also  shown  are  a  few 
points  calculated  from  (5.2).  To  plotting  accuracy  they  are  the  same  as  the 
solution  to  (2.2)  even  though  the  value  of  e  is  not  very  small.  Thus 
diffusion  can  be  included  analytically,  giving  results  essentially  the  same  as 
(2.2).  The  front,  or  shear  discontinuity,  is  at  S  *  0  and  from  (5.2)  its 

thickness  is  0(o^2)  or  O(E^). 


VI.  CONCLUSIONS 

The  Wedemeyer  model  for  spin-up  of  a  fluid  contained  in  a  cylinder  which 
is  impulsively  rotated  about  its  axis  has  been  implemented,  and  numerical 
solutions  to  the  spin-up  equation  have  been  obtained.  The  approximation  to  an 
impulsive  start  was  discussed.  A  critique  of  the  Fkman  compatibility  con¬ 
dition  based  on  the  Rogers  and  Lance  boundary  layer  mass  flow  relation  was 
given;  there  are  compelling  reasons  for  not  using  it. 
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The  finite  difference  method  of  solving  the  spin-up  equation  is  presented 
in  detail  and  the  results  of  error  studies  given.  Error  due  to  a  disconti¬ 
nuity  in  boundary  conditions,  caused  cy  the  impulsive  start,  were  avoided  hy 
developing  a  local  solution. 

A  number  of  results  were  presented  t.o  illustrate  the  main  features  of  the 
theory  and  calculations.  Comparison  of  the  results  with  those  from  finite 
difference  solutions  to  the  Navier-Stokes  equations  validate  the  spin-up 
theory  and  numerical  method.  The  effects  of  including  diffusion  art.-  ade¬ 
quately  shown  by  various  examples  hut  an  analytical  approach  is  also  derived. 
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Figure  1.  Notation  and  Coordinate  System  for 
the  Spinning  Cylinder. 
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Normalized  In-Bore  Spin  History  for  an  M687 
Shell  Launched  From  the  M109  Howitzer  at 
Charge  3,  s?  =  92.3  Rev/s. 
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Figure  3.  Normalized  Radial  Mass  Flux  in  the  Boundary  Layer  vs  Ratio  of 
Outer  Flow  Spin  to  Disk  Spin,  Reference  13.  Also  U  vs  V/r  when 
the  RL  Curve  is  Used  to  Define  an  EC. 
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Figure  4.  Grid  Scheme  in  (r,t)  Plane. 
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Figure  5.  V/r  vs  r  for  Re  =  4974,  c/a  =  3.30  at  t  =  20  and 
50  Using  Local  Solution  or  Assuming  V(l,0)-1. 
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Figure  6.  Vorticity  c  vs  r  for  Re  =  4974,  c/a  =  3.30  at  t  =  100 
and  600  Using  Local  Solution  or  Assuming  V(1,0)  =  1. 
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Figure  9.  V  vs  r  for  Case  2  With  c/a  =  5.200,  Re  =  2  x  10  at 
Three  Values  of  t;  lst'  2700,  Turbulent  Ekman  Layers 
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Figure  10.  V  vs  r  for  Re  =  3  *  10  at  Two  Values  of  c/a,  t 
With  Laminar  or  Turbulent  Ekman  Layers. 


Figure  11.  V  vs  r  With  Diffusion,  (2.2),  and  Without  Diffusion, 

(2.4),  for  Re  =  6.08  *  105,  c/a  =  2.679  at  Two  Values 
of  t.  Laminar  Ekman  Layers. 


Figure  12.  V  vs  r  With  Diffusion,  (2.2),  and  Without  Diffusion, 

(2.4),  for  Re  =  2  x  io6,  c/a  =  5.200  at  Two  Values  of 
kft;  kj.  =  l/t$t;  =  1/2700,  Turbulent  Ekman  Layers. 
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Figure  15.  V  vs  r  for  e  =  1/50  at  t'  =  0.80  From  (2.2)  and  (2.4). 

The  Venezian  Approximation  is  the  First  Term  of  (5.2). 
The  Points  are  From  Both  Terms  of  (5.2). 
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APPENDIX  A:  COEFFICIENTS  IN  EQ.  (3.7) 

=  -k jj  [j  Ar  -  (1/2)  (Vj  k  +  Vj  k+1)]  —  laminar  F.kman  layer 

=  -kt  [j  at  -  (1/2)  (Vj  k  +  vj ,k+l)j^J  "■*  turlju1ent  Ekman  layer 

=  ( At/ Ar)u  ,  D  =  [At/(Ar)2]/Re 

aN  “  0 

-(A./4)  +  [(4(j-l)}  ^  -  (l/2)]0  for  j  *  1,N 

1  +  (Aj/2)  +  (30/2) 

1 

1  +  [Aj/(2j)]  +  D  for  j  *  1,N 

Uj/4)  -  (3D/4) 

0 


(Aj/4)  -  [  {4(  j4-!)  )-1  +  (1/2)30 


for  j  *  1,N 


dl  "  vl,k  •  <V4)  <v2,k  +  2vl,k}  *  (D/4)  (3v2,k  "  6vl,k> 

dN  -  1 

dj  ’  Vj-l,k  {(Aj/4)  +  (D/4)  tMJ-1)"1])  +  vj|k  U-n-(Aj/[2j])} 

+  v  uk  {-( Aj / 4)  +  (0/4)  [2+( j+1)"1])  for  j  *  1,N 
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LI  SI  or  SYMBOLS 


cross-sectional  radius  of  cylinder 
half-height  of  cylinder 


2 

s  v/ a  c ' ,  Lkman  number 

abbreviation  for  "Lkman  compatibility  condition" 

index  indicating  r-coordinate  of  nodal  point  in  finite 
difference  solution  (r  =  j  Ar) 

index  indicating  time  in  finite-difference  calculation 
(t  =  k  At) 


s  <  (a/c)  Re"1/2  (see  (2.8)) 


'  0.035  (a/c.)  Re 


-1/5 


(see  (2.4)) 


-1/2 

radial  mass  flux  in  boundary  layer/r  Re 

number  of  subintervals  in  r  in  finite-difference  solution 

pressure/ ( pa2id2)  in  Wedemeyer  model  (see  (2,3)) 

zeroth  order  approximation,  in  ( 1 / 1 s ) ,  to  P*  (see  (2.6)) 

coefficient  of  first,  order  approximation  to  P*  in  (2.6) 

pressure/ ( pa" sV')  in  solution  to  Navier-Stokes  equations 
(see  (2.1)) 

radial  coord in ate/ a 

=  (1  -  r)/t,1/?- 

~  a2tVv,  Reynolds  number 

abbreviation  for  "Rogers-larice" 

(y/£)1/z 
t  i  me  x  ol 


st/ts 

time 

time  required  for  spin  of  cylinder  to  reach  *1;  i.e., 
~1  (tc)  =  £1 
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characteristic  time  of  flight  of  projectile 
’-1/2 

=  (d  si/dt)  ,  acceleration  or  Impulse  time 
1/2 

=  (2  c/a)  Re  /n,  characteristic  spin-up  time  for  laminar 
Ekman  layer 

1/5 

=  (28,6  c/a)  ke  /il,  characteristic  spin-up  time  for  turbulent 
Ekman  layer 

=  ( a/c ):?/(En),  time  scale  for  viscous  diffusion  of  vorticity 
=  2n/n,  time  for  one  revolution 


U,  V,  W 


U  ,  V 
w  w 


ur  Vlf  Wj 


V  V1 


radial,  azimuthal,  axial  velocity  components  x  l/(an)  of 
Wedemeyer  model  spin-up  flow  with  diffusion  (see  (2,2)  and 
(?..3)) 

radial  and  azimuthal  velocity  components  x  1/ (an)  of 
Wedemeyer  model  spin-up  flow  without  diffusion  (see  (2,4)) 

coefficients  in  first  order  approximations,  in  1/t  ,  to 
U*,  V*.  W*  (see  (2.6))  s 


U*,  V*,  W*  radial,  azimuthal,  axial  velocity  components  x  l/(an)  of 
Navier-Stokes  flow  (see  (2.1)) 


zeroth  approximation  to  V*  (see  (2.6)) 

functions  of  R  in  early  time  solution  in  Chapter  IIIC 

«  (rVV-  1)/(SV-  1) 

i 

axial  coordinate  (z  =  0  at  cylinder  midplane) 


=  rV,  circulation 

=  1/N,  r-ifiterval  size  in  finite-difference  solution 
t-interval  size  in  finite-difference  solution 


l/(k4  Re) 

=  (rV)p/(2r),  nondimensional  vorticit.y 
azimuthal  angle 

constant  in  expression  for  radial  velocity  with  laminar 
Ekman  layer  (see  (2.8)) 
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v  kinematic  viscosity  of  fluid 

p  density  of  fluid 

Si  angular  velocity  of  spinning  cylinde 

Fi  (t)  time  history  of  projectile  spin 
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